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ABSTRACT 

We present the first detections of the 3 /im water ice and 3.4 /im amorphous 
hydrocarbon (HAC) absorption features in z ~ 2 ULIRGs. These are based on 
deep rest-frame 2-8/im Spitzer IRS spectra of 11 sources selected for their ap- 
preciable silicate absorption. The HAC-to-silicate ratio for our z ~ 2 sources is 
typically higher by a factor of 2-5 than that observed in the Milky Way. This 
HAC 'excess' suggests compact nuclei with steep temperature gradients as op- 
posed to predominantly host obscuration. Beside the above molecular absorption 
features, we detect the 3.3 /im PAH emission feature in one of our sources with 
three more individual spectra showing evidence for it. Stacking analysis suggests 
that water ice, hydrocarbons, and PAH are likely present in the bulk of this 
sample even when not individually detected. The most unexpected result of our 
study is the lack of clear detections of the 4.67 /im CO gas absorption feature. 
Only three of the sources show tentative signs of this feature and at significantly 
lower levels than has been observed in local ULIRGs. Overall, we find that the 
closest local analogs to our sources, in terms of 3-4/iin color, HAC-to-silicate 
and ice-to-silicate ratios, as well as low PAH equivalent widths are sources dom- 
inated by deeply obscured nuclei. Such sources form only a small fraction of 
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ULIRGs locally and are commonly believed to be dominated by buried AGN. 
Our sample suggests that, in absolute number, such buried AGN are at least an 
order of magnitude more common at z ~ 2 than today. The presence of PAH 
suggests that significant levels of star-formation are present even if the obscured 
AGN typically dominate the power budget. 

Subject headings: galaxies:high-redshift,galaxies:ISM, infrared:galaxies,quasars:absorption 
lines 



1. Introduction 



One of the most exciting recent discoveries in astronomy is the strong evolution and 
co-evolution of both the s tar-formation rate density, and black hole number density from 
z ~ 1-3 to today (see e.g. IWall et alJboOfit Irlopkins fc Beacoml 120061 : Irlopkins et allhoosl ). 
Particularly interesting are composite or transition objects where quasa r s co-exist with pow - 
erful starbursts (see e.g. ICanalizo fc Stockton! 1200 ll ; ICoppin et al.ll2008l ; IWalter et al.ll2009l ). 
Such objects are likely to be much more common at z~2, the peak epoch of both the 
quasar number d ensity and the star -formation rate density. In particular, recent Spitzer 
Space Telescope ( jWerner et al.l I20Q4J ) results suggest that the ULIRG n umber density wa s 
more than 200-fold greater at z~2 than it is today (for a review see ISoifer et al.ll2008l ). 
This strong evolution of the luminosity function means that many more extremely lumi- 
nous sources (Lir, ~ 10 13 L ) have been discovered than are observed locally. By necessity 
therefore these extreme sources are typically compared with local ULIRGs which are in fact 
much less luminous (typically L IR ~ 10 12 - 10 124 L ). Some of these z ~ 2 sources derive their 
luminosities from powerful obscured quasars, but o thers appear to be scaled - up versions of 



the local largely starburst-dominate d ULIRGs (e.g. ISajina et al.ll2007al . 120081 ; iPolletta et al. 



20081 ; iHernan-Caballero et al.l 120091 ) . Due to the dearth of such sources nearby the physical 
conditions such as the dust composition at these highest luminosities is poorly understood. 

One approach is to determine whether or not the detectability and relative strength of 
absorption and emission features arising from different dust components are consistent with 
those of local ULIRGs in order to test how similar are their dust properties. The mid-IR 
regime is well suited to such a study as it contains a wealth of molecular absorption features. 
By far the strongest is the 9.7/im silicate absorption feature (Si-0 stretching mode). This 
is thought to sample the large silicate grains, which are frequently coated in various ices 
including water. These ices can be seen in the H-O-H bending mode (3.05 /mi) and the 
O-H stretching mode (6.15 /mi). In more extreme environments such as observed toward the 
Galactic Center, the degree of hydrogenation is such that prominent hydrocarbon (HAC) 
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features are observed (e.g. iPendleton et al.lll994t IWhittetll2003l ). The strongest of these is 
the 3.4 jum C-H stretching mode, but also includes the 6.85 and 7.25 /im C-H deformation 
modes. Studies of these features in external galaxies and in part icular ULIRGs are still quite 
challenging. The first such study was that of ISpoon et al.l (120021 ) who, out of a sample of 100 
sources with JSO-SWS spectra, detected the 6 /im water ice absorption in 18 galaxies (includ- 
ing most of the ULIRGs). Indeed ULIRGs show the de epest silicate absorption features (e.g. 
Hao et al.ll2007l ; ISpoon et al.ll2007l ; ISirocky et al.ll2008l ). and hence are the best candidates to 



look for the other mid-IR molecular a 



spectra of 37 ULIRGs (llmanishi et al 



jsorptio n features. The Subaru study of the 3-4/zm 
2006a) supplemented by the ^5-30 Aim Spitzer In - 



frared Spectrograph (IRS; iHouck et al.l 120041 ) study of local ULIRG s (llmanishi et al.l 120071) . 
and the VLT study of the 3-4/im spectra of 11 more ULIRGs by iRisaliti et al. I J2006h all 
show molecular absorption fe atures. Ove r the 3 - 5 /im regime, the VLT ISAAC spectra of 



5 ULIRGs were presented in ISani et al.l (120081 ). Most rec ently AKARI spectr a of 45 lo- 



cal ULIRGs provide the first continuous 2.5- 5 /im spectra (llmanishi et al. 



20081 ) of a large 



number of ULIRGs and serve as an excellent local basis for comparison with our high- z 
spectra. Lastly, CO gas has a prominent 4.67 /xm absorption featur e , whi c h has been de- 



tected in several deep ly obscured local ULIRGs (ISpoon et al.l 120041 . 120051 ; llmanishi et al. 
20081 : ISani et al.l 120081 ) . An archival study of ISO spectra of Seyfert galaxies showed no sign 
of this fe ature, presumabl y due to different obscuration geometries, i.e. dust torus vs. buried 
nucleus ( jLutz et al.l 120041 ) . All these absorptio n features arise in slightly different conditions 
and hence sample different obscuring layers (ISpoon et al.l 120041 ) . However, how these fea- 
tures arise and their relative strengths are still not well understood even locally. Only a few 
dust models th at address the relative strength of the differ ent absorption features have been 
published (e. g. lOssenkopf fc Hennineill994l ; IZinoveval 120051 ). which along with laboratory ice 
spectra (e.g. iGerakines et al.l 119961 ) help in interpreting the results. Most of what we know 
about interstella r ices comes from observations of embedded protostars (for a review see 
Gibb et al.l 120041 ) which in may also be analogs to our most highly obscured z ~ 2 sources 
than local ULIRGs. 

In this paper, we present unprecedented signal-to-noise rest-frame ~ 2 - 8 /im spectra of 
z~2 ULIRGs. A spectrum of a z~2.56 galaxy with similar signal-to-noise is presented in 
Rigby et al.l (120081 ). but starting at restframe ~ 5.5 /im. Spectra extending to the rest-frame 
~ 2 /im have several advantages which will shed light on the physical conditions of these 
distant systems: 1) the continuum arises from hot dust near the power source, 2) this regime 
includes the 3.3 /im and 6.2 /im PAH features - both isolated and less severely affected by 
the mid-IR continuum model than the longer wavelength PAH features, 3) dust absorption 
features including water ice at 3 /im and 6 /im, several hydrocarbon (HAC) features, and 4) 
the warm CO gas 4.67/im absorption feature. In this paper we concentrate on the various 
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absorption and emission features in these spectra and compare our results with the avail- 
able local ULIRG data. Throughout t his paper we assum e a flat Universe with f2 M = 0.27, 
Q A = 0.73, and H = 71km s" 1 Mpc" 1 dSpergel et al1l2003h . 



2. Sample 
2.1. Background 



Our group has been engaged in a multifaceted study of the properties of high-z 
ULIRGs, including low-resolution IRS spectra, HST imaging, X-ray imaging, and Keck spec- 
troscopy. We have ~200 IRS spectra of sources with F 2 4>0.8mJy, selected in the Spitzer 
First Look Surve y fielc0, of which ~ 70 sources are at 1.5 < z < 2.5, and likely L m ~ 10 13 L 
(jYan et al.l 120071 ). The anal ysis of the first 50 sources (GOl pro gram #3768) of this sample 
has already been completed (jYan et al.ll2007l ; ISajina et al.ll2007al ). The remaining 150 sources 
(G02 program #20629) give a largely uniform, flux-limited sample (Dasyra et al. in prep.). 
Both the GOl and G02 samples have a limit of F 2 4 > 0.8 mJy, but our GOl sources have two 
additional selection criteria: uF u (24 fim)/uF u (0. 64 /mi) > 1.0 and ^(24)/*^ (8) >0.5. The 
later selection cut means that the GOl sources tend to be fainter at 8 /im (rest-frame 2 /an) 
than the G02 sources. 



2.2. Target selection 



We adopt three selection criteria for our follow-up program: 1) redshift, 2) brightness, 
and 3) silicate feature depth. We adopt the redshift range 1.6 < ,2 < 2.5. These boundaries 
are determined by practical reasons: beyond z>2.5, IRS spectra poorly sample the 9.7 /im 
silicate feature, while sources below z~1.6 require an observationally-expensive additional 
IRS module (SL2). For the brightness, we simply require an 8 /im detection (>60/iJy) or 
else the SL spectra would be too difficult to obtain. The last criterion was most difficult to 
define. As this is a largely exploratory program, we did not want to restrict ou rselves to 
only the most extremely obscured sources (see e.g. category 3; ISpoon et al.l 120071 ) . but also 
needed the sources to have relatively strong silicate absorption features since, as expected, 
only sources with sil icate absorption have been observed to have ice absorption features 
(jlmanishi et al.ll2007l ). We selected, r 9 7 > 2.0 (based on the Galactic Center extinction curve 
used in Sajina et al. 2007a) as a rather generous lower limit. In terms of the observed 9.7 /im 



1 http: / / ssc.spitzer.caltech.edu /fls/ 
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depth (i.e. ln(F cont /F b s ) as is more commonly used, this translates to Tsi = 19.7 > 1.45. For 
the sake of comparison with previous work, throughout this paper we use observed optical 
depths (see also § !4.3p . Out of 52 sources in the GOl sample, 6 sources satisfy the above 
conditions. Of the G02 sample, 5 sources satisfy the conditions. However, as we used very 
preliminary data reduction of this sample, when doing the initial sample selection, one of the 
sources that meets the criteria (MIPS509) was not included, and another (MIPS22633) was 
observed although its silicate feature depth is too low compared with the rest (r 9 .7 ~ 1.1). In 
total, our sample consists of 11 sources. 



2.3. Redshift accuracy 



Only 3/11 of our sources have optical/near-IR redshifts (see Yan et al. 2007 and Sajina 
et al. 2008). The rest are typically based on the 10 /im silicate absorption feature and 
hence are believed to have an uncertainty, Az^0.2 (see Yan et al. 2007). Hence, for the 
sources without accurate redshifts, we look for example for absorption features within the 
range zirs ± Az. This constraint makes the reliability of single marginally-detected features 
doubtful. In cases of unambiguous and/or multiple absorption feature detections, these 
can be used to improve our redshift estimates. As the rest-frame peak of the molecular 
absorption features can also v ary - for example the 3 Aim ice f eature can peak anywhere 
in the range ~3.05- 3.1 /im (IChiar et al.l 120021 ; IZinoveval 120051 ). an uncertainty of up to 
Az ~ 0.05 remains. 



3. IRS observations 



In terms of both the molecular absorption features and hot dust continuum, we are 
interested in obtaining deeper spectra in the rest-frame ~2-8/iin spectral range. At the 
redshifts of our targets this means using the Spitzer IRS's SL1 and LL2 modules. One of our 
targets, MIPS22432 is very bright (F 8 = 0.5mJy), and was already observed in SL1 as part 
of our G02 program. Its signal-to-noise is sufficient for our purposes (see below) and hence 
does not need to be re-observed in SL1. Both MIPS22432 and MIPS15880 have sufficiently 
high signal-to-noise and do not need new LL2 observations. Lastly, MIPS8392, was already 
observed in SL1 as part of the GOl program #15 (P.I. Houck) and hence its SL1 observation 
was removed from our observations list. Unfortunately, the depth of the previous observation 
is 3 x lower than we needed here, but this is still considered a duplicate observation (the 
limit is 4 x). The depth of the earlier observation is insufficient for our purposes hence it 
can serve only to constrain the continuum levels, but is not useful in looking for absorption 
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features. 

Table [T] gives a summary of the observations covering three Spitzer observing cycles 
(G01,G02, and G04). The full GOl sample and observations are discusse d in detail in 



Yan et al.l (120071 ). The full G02 sample and observations are discussed in iDasyra et al. 
( 120091 ). The G04 observations (PID$40793) are presented here for the first time. Combin- 
ing all these programs, each source has typically ~2.5hrs on source of IRS SL1 and LL2 
observations, although brighter sources have less exposure time as can be seen in Table [TJ 
The SL1 observations were mostly in mapping mode in order to improve the signal-to-noise 
ratio in these very faint sources. For our sources, we adopted the IRAC and/or VLA po- 
sitions which are the most accurate available and sufficiently accurate for the SL1 module 
which has the narrowest slit width (3.6" compared with 10.5" for LL2). 



3.1. Data reduction 



We use IrsLow, a data reduction code for IRS low resolution observations custom-written 
by Dario Fadda (see Fadda et al. 2010 in prep, for further details). The data reduction 
steps include bad pixel removal, background subtraction, co-adding different positions (such 
as the 2 nod positions in stare mode, or the 4 positions we typically used in map mode), 
and finally extracting the spectra using optimal extraction. Bad pixels are determined as 
outliers across the available frames. For the SL data for example, fainter sources have 40 
frames per source, but the bright sources have only 6 frames per source. Therefore this bad 
pixel flagging is less efficient for our bright sources with very few individual frames leading to 
somewhat noisier spectra than the nominal difference in exposure time. The code is iterative 
in that after the first background subtraction, any sources present are masked, and a second 
background subtraction is done. In addition, we test the reliability of any features we see, by 
displaying all points that go into the spectral extraction and manually clipping any obvious 
outliers (on top of the bad pixel flagging of the code) that might bias our results. A few 
spurious spikes are removed in this way, but this is mostly a concern for narrow features of 
interest such as the 3.3 jum PAH feature. In all cases, where we believe that this feature is or 
might be present we have used this option of the code to check for any possible outliers. For 
the broad water ice or HAC absorption features, the quality of the background subtraction 
is more of a concern than any remaining outlier pixels. Below we describe our background 
subtraction quality test. 

As a double check on our per pixel errors (and the background subtraction) we also 
extract sky spectra from a position in the data image which does not include a source. The 
same PSF is used for the off-source spectra extraction as the source spectra. These sky 
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spectra allow us to ascertain that the mean background level is zero and that no residual 
gradients are present. The mean rms values from these sky spectra are found to be compara- 
ble with the typical per pixel errors of the extracted spectra (but ~ 30% worse in the case of 
the two SL stare mode observations due to the less efficient bad pixel flagging). For the SL1 
spectra the rms is typically ~ 20 - 50// Jy (with the lower values typical of the fainter sources 
and higher values for the brighter sources) with continuum SNRs of ~4-8. For the LL2 
spectra the rms is typically ~ 60 //Jy with continuum SNRs ranging from ~ 3 for MIPS22530 
to ~ 25 for MIPS22303, with typical values of ~ 10. 

We tested IrsLow's spectral extraction, by running SPICE after the background sub- 
traction and obtained fully consistent results. We find that our old reduction of the LL2 data 
for the GOl sources (jYan et al.l 120071 ) is frequently lower than the new data by up to 30% 



(somewhat worse for MIPS464). However, when re- reduced the same way as the new data 
(see above) the two sets are consistent. We attribute the difference to improved background 
subtraction and calibration. For consistency we have re-reduced all available data (including 
LL1, LL2 and SL1) before co-adding into the final spectra. 



4. Results 

4.1. Individual spectra 

Figured] shows the eleven spectra of our sources derived by co-adding all available data 
(see Table [TJ. For the LL2 region (A Q b s ~ 14-21 //m), the spectra represent the weighted 
mean of the old and new data. The panels below show closer views of the spectra around the 
principal water ice absorption features at 3 //m and 6 //m. We detect the continuum down to 
~2.5 //m allowing us to look for the various molecular absorption features expected in this 
regime. Indeed, we see a number of interesting features in the spectra, the vast majority of 
which can be identified with the known molecular features including water ice, PAH, HAC, 
and CO. Figured] also shows the SL1 and LL2 sky spectra for these sources (i.e. spectra 
taken off-source, but with the same PSF). As expected, these sky spectra have a median of 
essentially zero, but do show some 'bumps' and 'dips' which in addition to the background 
noise, likely also include some low-level rogue pixels that are not sufficiently strong outliers 
to have been removed. Such spurious features, however, could influence our source spectra 
as well. Essentially there are three considerations for believing the reality of a given feature: 
1) whether or not it is stronger than the amplitude of the features in the sky spectra, 2) 
whether or not it has consistent redshift, and 3) whether or not it has roughly the expected 
shape (i.e. width). We address these issues in § 14.31 Before hand however, we examine briefly 
the hot dust continuum shape of the sources. 
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Fig. 1. — The combined ~ 2- 12 /mi final spectra for our sample. The yellow bands denote 
the border regions where the different orders are stitched and hence any 'features' that 
appear in those bands are likely artefacts. A smoothing of 5 pixels is applied to all spectra. 
The pink shaded region denotes the 1 cx/pixel uncertainty in the spectra. The blue diamonds 
are the IRAC8 /im and MIPS24 /xm broadband fluxes respectively. The two panels below 
are close-ups of the regions around 3 /im and 6 /im, where the errorbar denotes the mean 1 a 
error. A sky spectrum is also shown (blue line). For easier comparison, this is often offset 
(F s k y +const.), where the constant in given in blue just above the sky spectrum. 
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Fig. 1.- 



continued. 
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Fig. 2. — Examples of continuum fitting in the 3-4/zm regime. Here we show all sources 
showing evidence of features in this regime. The crosses represent the actual spectra, while 
the solid red line is the spectra smoothed by a factor of 5. The dashed line is the best-fit 
continuum. 



4.2. Continuum slopes 

We can see in Figured] that there is considerable variety in the continuum shapes for 
our sources. This likely is indicative of differences in obscuration geometry. To address this 
fully requires comparison with radiative transfer power source+obscuration models which is 
beyond the scope of this paper. Here we merely wish to quantify the basic characteristics of 
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Fig. 3. — The 3-4 jU.ni spectral in d ex (d efined as f v (xv~ a ) for our sample and the AKARI 
ULIRG sample of llmanishi et al.l (120081 ). Many of our sources show steeper slopes than 
observed locally. 



these spectra such as spectral slope and monochromatic luminosity and compare with those 
of local ULIRGs. 



We adopt the standard measure of spectral slope which is a power law spectral index 
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(i.e. F v (x.v a or A Q ). We define continuum regions free of major features (2.6 < A < 2.8 /im 
and 3.8 < A < 4.2 /mi). We fit the above model to the spectra that fall within the continuum 
regions (see Figure^ . For MIPS464, MIPS15880, and MIPS22432, their lower redshift means 
that the spectra do not reach rest-frame 2.8 /mi. For these sources we adopt A < 3.2 /mi as 
the continuum region. This assumes there is no water ice absorption, which is not observed 
for MIPS464 and MIPS16152 and is not possible to constrain for MIPS15880 or MIPS22432. 
We find fairly red spectra in this regime with slopes in the range a ~2-4. 

To compare our values with those of local ULIRGs wi thout the bias o f differ ent contin- 
uum definitions, we take the ULIRG AKARI spectra from llmanishi et al.l (120081 ) and apply 
the same continuum fitting to them. Figure[3] shows the distribution of the 3-4 /im slope for 
our sample compared with local ULIRGs. We find t hat the bulk of our sample has steeper 
slopes than the local ULIRGs. iRisaliti et al.l (120061 ) argue that a red 3 -4 /mi slope is an 
indicator of an AGN dominant in the mid-IR (although no t necessarily bolome trically domi- 
nant). As defined, a is related t o the slope V a s defi ned by IRisaliti et al.l (120061 ) via T = a-2. 
Using this we can translate the IRisaliti et al.l (120061 ) criterion to: a > 2 sources are predom- 
inantly AGN-dominated. This is satisfied by most of our sources, wh ich is not surprising 
as most are already believed to be AGN-dominated (jSajina et al.l 120081 ) . MIPS22530 is our 
strongest- PAH, starburst-dominated source and it has the flattest spectrum with a = 1.9. 

We follow the same procedure for the rest-frame ~ 5 - 7 /im spectra, where we define 
as continuum regions 5.0<A<5.5/im and 6.7 < A < 7.0 /im. The slopes there tend to be 
considerably flatter (except where strong-PAH contaminate this slope) as can already be seen 
visually in Figured] Lastly, we compute the monochromatic 3 /mi and 5.5 /im luminosities 
based on these fits. Table[2] gives the continuum parameters (luminosites and slopes) of our 



sources. 



4.3. Optical depth measurements: feature reliability 

In the previous section, we discussed the determination of the continuum levels for 
our sources. Figure[2] shows the best-fit continua for the sources with water ice and/or 
3.4 /tm HAC absorption (MIPS464, MIPS16059, MIPS16080, MIPS22303, and MIPS22633). 
This figure also shows that the absorption features are not driven by smoothing over a 
single outlier/spike. None of the sky spectra of these sources show broad features of such 
strength which also excludes them being background subtraction artefacts. After the con- 
tinuum fitting, the optical depth of the various features is then determined based on: 
r = ln(F cont /Fdata)- We should stress that here we use the convention of calling this ob- 
served feature strength an optical depth, even though the true optical depth is dependent 
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Fig. 4. — The 3.05 /zm water ice and 3.4 /xm HAC opacity for the sources where these are 
observed. We have applied a binning of 5 pixels to all spectra to increase the clarity and 
effective signal-to-noise {rms / V 'Nun) ■ The dashed and dotted lines are the ave rage profiles 



of the 3.05 /iin water ice and 3.4 /zm HAC features for Galactic center stars from lChiar et al. 



Chiar et al. 


2002; 


Dulev et al. 


2005) 



hence a wide range of specific shapes is expected (e.g. 
MIPS16080 and MIPS22633 show evidence of both ice and HAC features but the combined 
profile is inconsistent with expectations if these are the only contributors. The observed 
spectra could be expl ained by the presence of the PAH 3.3 jttm emission feature as well (e.g. 



Imanishi et al.ll2006al ). 
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on the exact radiative transfer solution. In particular, in the case o f the realistic clump y 
dust distribution the two do not translate in any straightforward way (INenkova et al.ll2008l ). 
As discussed previously, the reliability of weakly detected features is complicated by the 
frequently uncertain redshifts (except for MIPS16059, MIPS16080, and MIPS22530 which 
sources have Keck rest-frame optical spectra). For MIPS464, MIPS22303, and MIPS22633, 
the observed features are of a redshift consistent with the silicate feature and sufficiently 
strong in themselves to make the feature identification convincing. Table[3] gives a summary 
of the detected features as well as the updated redshifts. The uncertainties are estimated 
based on Fdata ± la. We ignore the uncertainty in the continuum level. 



4.4. Water ice and HAC absorption 

In § !4.3l we claim detections of water ice and/or HAC absorption for 5/11 of the sources in 
our sample (see Table [3]) • Another test of the reliability of these features is to compare their 
observed profiles with the average profiles of these features as determined for Galactic sources. 
In FigureHJwe show the same five sources as in Figure[2]but with the continuum subtracted, 



and th e average profiles of the 3.05 /iin water ice and 3.4 /im HAC features from lChiar et al. 



( 120021 ) overlaid. The continuum is the same as determined in the previous section, although 
it should be kept in mind that there is some uncertainty in the continuum-fit parameters, 
as the fit is done on fairly small sections of essentially feature-free spectra. We find that 
there is broad agreement although our spectra show noise about these profiles. The ~ 3.0 /im 
absorption feature in MIPS16059 is at somewhat lower wavelengths than expected for water 
ice and therefore we treat this as an upper limit only and not as a detection. Overall, 
the most significant deviation is observed for MIPS16080, MIPS22303 and MIPS22633 all 
of which show an excess at 3.3 /im in addition to the likely water ice and/or 3.4 /im HAC 
absorption. To test the possibility of this being due to the 3.3 /im PAH, in Figure O we 
compare the MIPS22633 and MIPS 16080 spectra with reasonably closely matching local 
ULIRGs where we degrade the resolution of the AKARI spectra to m atch that observed. 



Both local comparison sources show 3.3 /im PAH ( llmanishi et al.ll2008l ). although in IRAS 
19254-7245 it is much weaker than in IRAS14121-0126. We find that indeed the 3.3 /im 
excess can very well be due to PAH emission especially for MIPS22633. For MIPS16080 the 
comparison is less convincing and therefore we treat its 3.3 /im emission as an upper limit 
only. 

For MIPS22303 we find an excellent local analog in the JSO-SWS spectrum of the Galac- 
tic Center source GC Ring NE (see Figure[6]). This comparison suggests, that MIPS22303 
is well described by a combination water ice+3.3/im PAH without any clear need for HAC 
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Fig. 5. — The spectrum of MIPS22633 compared with the AKARI spectrum of IRAS14121- 
0126 (z — 0.151). The later has been degraded to the spectral resolution of the MIPS22633 
spectrum. The combination 3.3 /im PAH+3.05/zm water ice absorption for MIPS22633 is 
quite feasible. 

absorption (as shown in FigureH]). For this source therefore, we only have an upper limit on 
the HAC absorption. 

4.5. 4.67 /im CO absorption 

One of the goals of this program was to try to detect the warm CO gas absorption 
feature at 4.67 /im. Unfortunately, in several cases, rest-frame 4.67/tm falls right in the 
region of where we stitch together the SL1 and LL2 orders (yellow band in Figured]) making 
any features observed there unreliable. Sources with z > 2.1 (i.e. where this region is entirely 
contained in the LL2 spectra) or with ,2<1.8 (i.e. where this region is entirely contained 
in the SL1 spectra) are the only feasible ones to look at. There are a total of 7 sources 
that satisfy these criteria. Of these, three sources show signs of absorption in this regime: 
MIPS464, MIPS16059, and MIPS22432 (see Figured]). These features are fairly convincing 




-16- 



CL 
0) 

TD 

a 
u 

O 



-0.2 




-0.4 



-0.6 



2.6 



21 



3.0 3.2 3.4 3.6 

Rest wavelength [microns] 



3.1 



4.0 



Fig. 6. — The unsmoothed spectrum of MIPS22303 (green curve) compared with the ISO- 
SWS-01 spectrum of the north eastern part of the molecular ring around the Galactic Center 
(black curve). For MIPS22303, the combination water ice and 3.3/im PAH without HAC 
absorption seems quite likely. 



in themselves (compared with the background), but do they agree with the expectations 
for the CO feature? In Figured we compare their smoothed normalized flux spectra (i.e. 
data/ cont), with the sam e for the hyperluminous galaxy F00 183-71 11 and the IRAS00397- 
1312 (ISpoon et al.ll2005l ). While the nominal position of the line is at 4.67 /xm, the peak 
absorption in our sources (assuming it is CO) is typically blueward of that. As the CO gas 
absorption feature is composed of two broad branches (R and P), if one of these is stronger 
than the other it results in the obse rved shift. A stronge r R branch (shorter- A) absorption is 



also observed in IRAS08572+3915 (jlmanishi et al.ll200fih and IRAS00397-1312 (ISpoon et al. 



20051 ; llmanishil 120091 ) . Note that for all these galaxies, the redshift is well established on 
the basis of other features, without any additional shifts required to match the CO feature. 
For MIPS464 the redshift is based on the silicate absorption plus 3.4 /im HAC feature. For 
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Fig. 7. — The three sources with signs of CO absorption. The y-axis shows normalized 
flux (i.e. t =ln(normalized flux)). For comparison we overlay the same for 00183-7111 and 
00397-1312. The noisiest spectrum of our sources is that of MIPS22432 which is responsible 
for the large spikes observed. 

MIPS16059, the redshift is based on Keck spectra, but is also consistent with the silicate 
and 3.4 /im features. For MIPS22432, the only strong- PAH source of the three, the redshift 
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is based on the 6.2 /mi PAH feature. 

The observed strength of the CO feature for these three sources (the strongest in our 
sample) is r 46 ~0.2 (as discussed in §4.3 this is just the observed feature strength and does 
not translate to intrinsic opti cal depth) . This is somewhat weak compared with the CO 



feature strengths observed by ISani et al.1 (120081 ). which range from ~ 0.1 -2.2. The relative 



weakness of our features can also be seen in Figured 

An possible explanation for both the observed asymmetry and the overall weakness 
of the CO feature might be the effects of a starburst-powered continuum originating from 
material outside the CO absorption layer. MIPS222432 (a strong-PAH source) has a clear 
star burst- component and as we show in §4.8, a starburst component is likely present in all 
sources even if not individually observed. 



4.6. PAH features 



We look for the 3.3 /mi and 6.2 /mi PAH features which are not as confused with th e 
continuum model as for e.g. the strongest PAH feature at 7.7 /mi (jSajina et al.l l2007al ). 
and hence might reveal t he suspected presen ce of star- formation even in sources dominated 
by nuclear activity. In (jSajina et al.l l2007al ) we showed that the stacked spectra of weak- 
PAH (based on the shallower spectra) sources show P AH emission. Ther e are several cases 
of weak- PAH sources that have cold dust detections (jSajina et al.l 120081 ) . In addition two 
our sources, MIPS22432 and M IPS22530 were iden tified as strong PAH sources based on 
our previously defined criteria (jSajina et al.l l2007al ). MIPS22530 was also shown to be a 
st arburst- dominate d source based on its rest-frame optical spectrum and cold dust detection 
dSaiina et aDl2008h . 



MIPS22432 already had the necessary signal-to-noise ratio and hence no new data for 
it is presented here. We merely re-reduced the spectra and again see a clear 6.2 /mi PAH 
feature. Unfortunately, the 3.3 /mi feature falls right on the edge of its spectrum and hence 
we cannot test for its detectability. 

MIPS22530, on the other hand, shows an unambiguous 3.3 /mi PAH feature (see Fig- 
ured]). In addition, three more sources MIPS16080, MIPS22303 and MIPS22633 show some 
evidence for the 3.3 /mi PAH feature but somewhat confused with water ice and HAC ab- 
sorption (see Figure]!}. For MIPS22303 and MIPS22633, a potential 6.2 /mi feature falls 
in the unreliable wavelength range where we stich the different IRS orders (see Figured]). 
MIPS16080, whose 3.3 /mi PAH feature is too weak to be considered a detection, does show 
a tentative detection of the 6.2 /mi feature. The PAH emission of this source is made more 
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reliable due to its having an accurate optical spectroscopic redshift. In addition, MIPS16080 
has sta rburst-like optical l ine ratios, but AGN-line IR spectrum suggesting it is a composite 
source (jSajina et al.ll2008l ). 



We fit a simple local continuum+Lorentz profile model to each of the detected or 
hinted at 3.3 /mi PAH features. Figure[8] shows the results. Wherever it was possible, we 
did the same for the 6.2 /im feature. For MIPS16080 we find fe.2 — 2.5 x 10~ 15 erg/s/cm 2 
(EW~120nm). For MIPS22530 we find / 3 . 3 ~ 0.97 x l(r 15 erg/s/cm 2 (EW~550nm) and 
/ 6 . 2 ~3.89 x 10~ 15 erg/s/cm 2 (EW~370nm). These imply a 6.2/3.3 ratio of ~3-4 for 



these sources. This is ro ughly consistent with what is observed for local ULIRGs (e.g. 



Imanishi et all 120071 . 120081 ) although the observed ratio has a large scatter. Note however, 
that for MIPS16080, both the 3.3 and 6.2 /im PAH features are somewhat unreliable (see 
e.g. Figure[5]) and should really be considered upper limits. 



4.7. Comparison with local analogs 

Figured shows our sources overlaid with local analogs. From top to bottom we show 
sources in order of increasing silicate feature depth. We picked IRAS1925 4-7245, and 



IRAS 17044+6720 due to their being among the few sources in the AKRI sample of llmanishi et al 



(120081 ) which had comparably red 3-4 /im spectra. Here they are compared with MIPS16059 
and MIPS464 respectively, to show that selecting on the basis of red 3 -4 /an spectra alone, 
yields good agreement across the full observed range (~ 2—11 /im). MIPS22633 looks quite 
extreme and unusual, showing a very strong ~ 6 /tin absorption complex, while being the 
weakest silicate absorption source in this sample. However, we do find a local analog, 
IRAS23515-2917, which has an even more extreme 6 /im absorption coupled with a rela- 
tively shallow 9.7/im silicate absorption. Lastly our strongest silicate absorption source, 
MIPS8392 se ems to have a mid -IR SED remarkably like that of the hyperluminous galaxy 



F00183-7111 (ISpoon et al.ll2004j ). This galaxy is one of the very few local examples that have 
comparable luminosities to those of our sample, and therefore finding such good agreement 
is not surprising. 



4.8. Average spectrum 

In the previous sections we have shown that a number of features are detected in the 
individual sources spectra including the 3.4 /im HAC feature as well as the 3.05 /im water 
ice absorption feature. There are hints of some other features in the sources such as marked 
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Fig. 8. — The fits to the 3.3 /zm PAH emission for the four sources that show signs of it 
this feature. The data points show the spectra themselves, the solid black lines show the 
smoothed spectra. The best-fit continuum is given by the blue dashed line, whereas the total 
best-fit model is given by the thick red line. 

in Figured], but they are typically of low significance. These include in particular the 6 fim 
absorption complex (likely including the 6.15 jum water ice feature), the 6.85 and 7.25 /zm 
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wavelength [|0,m] 

Fig. 9. — From top to bottom, in order of increasing silicate feature depth: MIPS22633, 
MIPS16059, MIPS464, and MIPS8392. The red and blue curves show respectively the Spitzer 
IRS and AKARI (when available) spectra for local analogues. From top to bottom these 
are: IRAS23515-2917, IRAS19254-7245, IRAS 17044+6720, and F00183-7111. 

HAC features and the 4.67 /im CO feature discussed in § 14.51 To address the question of 
whether or not these are present in our sample, we stacked the spectra as shown in FigurefTUl 
Because there are too few sources to test the robustness of the stack by taking randomly half 
of it, we used a different technique which involves removing each existing pair of sources, 
or 28 possibilities. This implies that features that appear in each of the 28 stacks need 
to be present in the majority of sources in order not to disappear in at least some of the 
combinations. 
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Fig. 10. — The red curve shows the average spectrum of all sources excluding the strong- 
PAH sources MIPS22432, and MIPS22530 and MIPS8392 which is of too poor signal-to-noise 
below ~ 5 (Ltm. To assess the effect of individual sources, the grey curve show average spectra 
containing 6/8 of the above sources. In each of the curve a different pair of sources is excluded 
in order to show that some aspects of the stack are there independent of which sources are 
included or not (see text for details). The stacks are performed on un-binned data with all 
individual spectra normalized at 5.5 /im before co-adding. This point was chosen because it 
is in the center of our new spectra, it is pure continuum in all sources unaffected by edge 
effects and near the 6 /im water ice absorption which we want to look for. The only clearly 
outlying set of stacks are the four ones that exclude MIPS464. 

FigurefTTl shows that in the stacked spectra, we see that the 3.4 /im HAC feature is quite 
prominent, unsurprising given that we detect it in four of the individual spectra. Potentially 
there is also 3.05 /im water ice absorption and as well as potentially the 3.3 /tm PAH emission 
given some 3.3 /xm excess. The 3.0 /im feature is badly fit at the lower wavelength end for 
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Fig. 11. — The stacked 3-4/^m optical depth profile overlaid with the ice and HAC profiles 
as in FigureHl To obtain this, a continuum was fit to each stacked spectrum. As in FigureH] 
the red line is the average spectrum with the grey lines being stacked spectra with pairs of 
sources removed. 

most, but not all permutations. MIPS16059 is probably responsible for some fraction of 
these (see § 4.4), the rest is probably due to edge effects in sources without any clear features. 
Figure[10l however, suggests that the 6 /im absorption complex is there in most sources. The 
most prominent such feature is seen in MIPS22633 where the 6 /im absorption is clear in the 
individual spectrum as well. One thing to keep in mind is that in a number of the sources, 
the ~ 6 /im features in particular are not individually detected due to their falling inside the 
unreliable IRS modules overlap region. This uncertainty, as well as the likely 3.3 /im PAH 
make it difficult to estimate the average 3.0 /im water ice and 3.4 /zm HAC feature depths. 
FigurefTTl suggests that these are (t 3 . ) ~ 0.22 and (T3.4) ~ 0.30, but it should be born in mind 
that these are fairly uncertain. 
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The stacked spectra show signs of PAH emission superimposed on the strong mid-IR 
continua (recall that we have excluded the two individually strong-PAH sources in the sam- 
ple) . The ~ 8 /im profile suggests the presence of the 7.7 jum and 8.6 /im features. In addition, 
the 3.3//m feature (seen tentatively in the individual spectra for MIPS16080, MIPS22303 
and MIPS22633) appears clearly present. The stacks show a possible 6.2 /im feature as well; 
however, caution must be taken as for several of the sources this feature falls in the unreli- 
able edge region. To test the reliability of this 'detection' we also stacked only the sources 
with reliable spectra in this region (i.e. MIPS464, MIPS15880, MIPS16080, MIPS16152, 
MIPS22548). The 6.0 /im absorption and the 6.85 and 7.25 /im HAC features appear in this 
as well, although the 6.2 /im 'bump' is less pronounced than seen in FigurefTUl 

The 4.67 /im CO feature is hinted at but seems to be driven by the feature in MIPS464 
as when this source is not included we have the four outlier stacks. This re-inforces the 
conclusion that, apart from the three sources discussed before, the CO feature is either not 
present or is extremely weak (74.6 << 0.2) in the majority of our sample. 



4.9. Silicate feature depth reliability 



The silicate featu r e does not fall within our new de eper data and i s disc ussed in more 
detail in ISajina et al.l (j2007al ) for the GOl sources and iDasyra et al.l (120091 ) for the G02 
sources. However here we address the measurement and reliability of the silicate feature 
depths for our 11 targets specifically. This is necessary because in order to put our newly 
detected molecular absorption features in context with local sources we would like to look at 
the ice-to-silicate and HAC-to-silicate feature ratios (which we do in the Discussion). Many 
of our sources have poorly defined continuum beyond rest-frame ~ll/iin and also tend to 
have very low signal-to-noise inside the feature itself (see F igured]) . T hese features were 



originally fit with a power law plus extinction law model (see lSajina et al.ll2007al ). which as 
discussed earlier relates to the observed silicate feature depth via T 97rnode i=lAr 97j0bs . For 
some of the sources this is fairly well determined, but for others the fitted value has large 
uncertainty (either due to poorly determined continuum, or uncertainty due to the PAH 
features flanking the feature). For example, for MIPS464 we derived model 79.7 = 4.9 ±1.6 
which translates to observed 79.7^ = 3.4 ± 1.1. By contrast, MIPS15880 has considerably 
better defined feature, and we derive model r 9 7 = 3.9 ±0.4 which translates to observed 
79.7,06s = 2.8 ±0.3. In the case of the strong-PAH source, MIPS22530, the silicate feature is 
saturated and we can only determine that it has observed silicate absorption of 79.7 > 3.7. 



In order to further test the uncertainty associated with our silicate feature optical depth 
estimates, we also fit all of them using the same procedure as used for the local ULIRGs in 
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Spoon et al.l (120071 ). We find that for the majority of the sources, the two approaches give 



consistent results, the values obtained differ most for the sources with poor continuum and 
noisy features, but are within the fitted uncertainty (e.g. MIPS464 has r Si = 2.5 consistent 
with our r 9 7 = 3.4 ± 1.1). The most substantial difference comes for MIPS22530 where we 
obtain r$i = 1.7, considerably shallower than the value obtained from our model fitting. This 
is a well known uncertainty arising from the different treatment of the PAH features. Since 
MIPS22530 and MIPS22432 (the other strong-PAH source) do not have either HAC or ice 
detections, they do not feature in our optical depth ratio estimates (see Discussion) and 
hence this substantial uncertainty does not affect our results. 



Discussion 



5.1. PAH properties at z ~ 2 



Figure [12] shows the EW3.3 vs. EW6.2 for some representative local ULIRGs. We 
also show the derived values for MIPS22530 and MIPS 16080 (which is an upper limit only 
as neither feature is unambig uously detected by itse lf). We show the starburst-dominated 
sou rce limits as determi ned by llmanishi et al.l (j2006al ) based on optical spectral classification, 
and lArmus et al.l (120071 ) based on neon line ratios. On both counts, MIPS22530 satisfies the 
starburst-dominated ULIRG criteria, consistent with our prior knowledge of this source. 
MIPS16080 on the other hand is more likely a composite source, as expected. 

Four of the sources showed signs of the 3.3 fim PAH feature, with MIPS22530 having the 
most unambiguous detection. Based on the fits presented in §4.6, we can estimate 3.3 /im 
PAH log luminosities of 9.6, 9.6, 9.8 and 10.1 for MIPS16080, MIPS22303, MIPS22530, and 
MIPS22633. The t otal 3 - 1000 /im infrared luminosities for the first three sources are given in 
Sajina et al.l (120081 ) and for MIPS22633 we have log(L IR /L ) = 12.9 (Sajina et al. in prep.). 
For our most reliable 3.3 Aim PAH source (MIPS22530) we have L 3 3 /L T r = 8.9 x 10" 4 . This 
is consistent with the L 3 . 3 /Lir = 8.5 x 10~ 4 reported by lSiana et al.l (120091 ) for a lensed z ~ 3 
Lyman Break galaxy, the only other high- z measurement of the 3.3 /mi PAH. T his ratio is 



somewhat larger than observed for local ULIRGs (~ 10~ 4 ; llmanishi et al.ll2006al ). However, 
this is believed to be an aperture effect, as this ratio is observed to increa se with distance 
with the more distant ULIRGs having a ratio of ~ 10~ 3 (IMouri et al.lll990l ). 



For our other sources we find L 3 . 3 /L 1R =6.6 x 10~ 4 (MIPS16080), 4.3 x 10" 4 (MIPS22303) 
and 1.6 x 10 -3 (MIPS2 2633). If we take a ratio of ~ 10~ 3 to be indicative of starburst- 
dominated systems (see ISiana et al.ll2009l ). then the above can give us an alternative look 
into what fraction of the Lir is due to starburst activity and also what SFRs are feasible. 
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Fig. 12. — The equivalent widths of the 3.3 /zm and 6.2 /im PAH features in local ULIRGs 
compared with MIPS22530 and MIPS16080 (red dots). The local ULIRG data as well 
as the criter i a for starburst-dominated sources are based on llmanishi et al.l (l2006al ) and 
Armus et all J2007h . 



Accordingly, MIPS 16080 likely has equal AGN and starburst contributions to the bolometric 
luminosity, whereas for MIPS22633, the starburst may well be dominant. In MIPS22303, 
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the AGN is likely the dominant component. These are consistent with our earlier results 
( ISajina et al.ll2008l ). Note however, that using the above translation between L 3 3 and Ljr 
implies starburst-only luminosities of >10 12 L Q even for MIPS16080 and MIPS22303, and 
consequently SFRs of ~100M o /yr. However, as discussed in §4.6, these features being 
of low signal-to-noise ratios and somewhat confused with the ice and HAC absorption fea- 
tures, are not particularly reliable. The stacked spectrum shows that PAH emission is likely 
present in the majority of these sources, even excluding the two strong-PAH sources. The 
fact that star-formation is present is further supported by the fact that although most of 
these sources are not individually detected in the far-IR, their stacked 1.2mm (rest-frame 
~400/im) flux is strongly detected at ~0.5 mJy which is comparable to fainter (S , 850 ~2- 
3mJy) sub-mm galaxies (ISajina et al.ll2008l ). These numbers suggests that SFRs on the 
order of 100 M Q /yr are feasible for the bulk of the weak- PAH sample. Although in itself 
significant, star-formation at this level is insufficient to dominate the bolometric luminosities 
of these sources. 



5.2. Ice and HAC detectability compared with local ULIRGs 



Local ULIRGs show a high level of detec tability of water ice absorption features. For 
example the 46 sources in the AKARI sample (jlmanishi et al.ll2008l ) were selected out of the 
IRAS 1 Jy sample merely on the basis of visibility to AKARI and are therefore unbiased with 
respect to the local IRAS ULIRG population as a whole. The 3 jum feature was detected 
in 30/46 sources (65%), whereas the 3.4 /xm HAC feature was only observed in 4 of the 
sources (9%). The presence of the HAC feature has been suggests as an indicator of AGN- 
dominance in the mid-IR (IRisaliti et al.l 120061 ; llmanishi et al.l l2006al ) and hence this low 
fraction is consistent with local ULIRGs being predominantly starbursts. To determine the 
appropriate detectability for our sample we need to consider first that MIPS8392 has too 
low a signal to noise in this regime, and three other sources are of z<1.8 which makes this 
feature too near the edge for reliable detections. This means that 3/7 (43%) of our sample 
show the 3 /im water ice absorption feature, this is roughly comparable to the local ULIRGs 
given the small number of sources in our sample. On the other hand (5/11) of the sources 
(45%) show the 3.4 /im HAC absorption feature, which fraction is significantly higher than 
observed for local ULIRGs (even accounting for the Poisson error in such a small sample). 
However, the 3.4 /im HAC feature can be practically difficult to detect in sources with strong 
3.0 /im ice and/or 3.3 /im PAH emissio n. This is one of th e explanation for the low HAC 
detectability among the local ULIRGs (jlmanishi et al.ll2008l ). 



We have already shown that the best analogs to our sources are the sources with red 
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3-4 um spectr a rath er than the more typical local ULIRGs. There are five sources in the 
Imanishi et al.l (120081 ) AKARI sample that meet our criteria (i.e. ct3_4>2). Of these, 3 



show the HAC feature, and a different set of 3 show the 3.0 /im ice feature. This is roughly 
in better agreement with the observed detection rates for our sample (but with very poor 
statistics!). 



5.3. HAC-to-silicate ratio 



FigurefTBI shows the HAC/silicate ratio vs. the silicate feature strength for our sources 
compared with local ULIRGs. To construct this figure , we use the T3.4 measurements for local 



ULIRGs from several studies ( Imanishi et al.l l2006al ; iRisaliti et al. 



2006: 



Sani et al. 



2008 



Armus et al.l J2007h and 



Imanis hi et al. 120081) combined with the silicate feature depths from 
Sirocky et al.l (120081 ) . We find that, while the HAC detectability of our sample is higher than 
the local ULIRG population, the HAC/silicate ratio appears consistent with that observed 
among local ULIRGs where the HAC feature is detected. 

For comparison we also show the data for embedded protostars from the review paper 
of Gibb et all ( 2004 ). Note that the HAC-to-silicate ratio observed in the MW is remarkably 
constant across a wide r ange of silicate feature depth. This is consistent with the diffuse 
ISM value of 0.05-0.06 (jPendleton et al.lll994l ). Moreover, a recent study of the spatial 
variation of the optical depth in the Circinus galaxy (the nearest AGN to the MW) found 
that altho ugh the absolute op tical depth changed dramatically, this ratio remained nearly 
constant ( Colling et al.l 120091 ) . An enhanced HAC-to-silicate ratio was first reported by 
Imanishil ( 120001 ) who studied a sample of 4 local obscured AGN nuclei including NGC1068 
(shown in Figure fTBI) . Their favored interpretation of this enhanced ratio comes from the 
following explanation: each of the two optical depth samples the column density between us 
and the sources. If the dust is well mixed and in a screen approximation is valid the ratio 
is expected to remain constant (reflecting the relative abundance of the silicate and carbon 
grains and their optical properties). However if a steep temperature gradient between the 
~300K dust at 9.7/im and ~1000K dust at 3.4/zm is present, since the hotter dust is 
spatially closer to the power source, the HAC feature arises from a larger column density 
than the silicate feature and he nce the enhanced rat io. This is what is commonly referred to 
as a 'buried nucleus' (see Fig.2 llmanishi et al.ll2007l ). Notice that this model addresses only 
the obscuration geometry: it is neither a statement about the nature of the power source 
(although commonly taken to be an AGN in these cases), nor about the absolute level of 
obscuration (the embedded protostars show comparable silicate feature depths). 



An alternative explanation might be that deeply obscured ULIRGs and our z ~ 2 sources 
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Fig. 13. — The ratio of the 3.4 /iin HAC optical depth to 9.7 fim silicate feature optical depth 
vs. the 9.7 /xm sil icate feature. For comparison we show the data for embedded protostars 



(Gibb etal 



f lGibb et al 



2004) and local ULI RGs (see text). We also indicate the location of SgrA* 



20041 ) and NGC1068 (jlmanishil |2000| ) . The orange hexagon shows the average 



value for the sample based on the stacked spectra (but see §4.8 for caution). 



simply have a higher proportion of hydrocarbons. For example, bombardment with H atoms 
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can hydrogenate amorphous carbon grai ns increasing the HAC-to-silicate ratio (see e.g. 



Grishko fc Duleyl |2000| ; iDuley et al.ll2005l ). Various other processes can be important, s.a. 



UV radiation processing of ice mantles can lead to refractory mat erials on the grain surfaces 
that lead to the 3.4 /im absorption feature (IGreenberg et al.lll995l ). consistent with the some- 
what weaker ice features. Ultimately, the degree of HAC absorption depends on the balance 
of processes that support hydrocarbon generation to those that suppress it (s.a. converting 
hydrocarbons to aromatics). It is not unreasonable to suppose that such processes might be 
more efficient in our sources than the typical Milky Way sources. The geometric explanation 
is simpler and therefore favored, however, a difference in composition cannot be excluded. 



5.4. Ice to silicate ratio 



Figured!] shows the ratio of the 3.0 /zm ice feature optical depth to the 9.7/im silicate 
dust optical depth. As before, we compare our sample with both local ULIRGs and the data 
for MW stars. We find that unlike the case for the HAC-to-silicate ratio which is remarkably 
similar across for all available MW sources, the ice-to-silicate ratio shows a wide spread of 
mor e than an order of magnitude . This ratio is st rongly dependent on the ice mantle thickness 



(.e.g. 



Ossenkopf fc Henninglll994l ; IZinovevall2005l ). but radiative transfer details following the 
arguments presented in the previous section likely dominate. When we consider only the 
sources with red 3-4/zm spectra, and where both ice and silicate optical depths are available, 
we find that IRAS08572+3915 has r 3 . /r 9 . 7 = 0.08, and IRAS 19254-7245 has r 3 .o/r 9 . 7 = 0.15. 
Although a much larger sampling is needed, this suggests (especially in conjunction with our 
own sample) that sources with redder 3-4/im slopes have lower observed ice- to- silicate 
ratios. 



5.5. The 4.67 /im CO feature 

The presence of CO gas in at least some of these systems is known independently through 
mm- wave spectroscopic observations (Yan et al., in prep.). In particular, MIPS 16080 and 
MIPS22530 both show significant masses of cold CO gas. Unfortunately both of these are 
among the sources, where the 4.67/zm feature falls in the unreliable parts of the co-added 
spectra. Moreover, the presence of CO does not automatically imply the presence of the 
4.67 /im CO feature which arises from the much warmer regions in the vicinity of the galactic 
nuclei. 

We would expect to see this feature in our sources as it is seen in some deeply obscured 
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Fig. 14. — The ratio of the water ice optical depth to silicate feature depth. Symbols are 
as in Figured^! The orange hexagon shows the average value for the sample based on the 
stacked spectra (see also §4.8). The dashed lines represent the approximate ratios expected 
for different ice mantle thicknesses (as indicated). Even when ice absorption is detected, our 
sources are consistent with fairly thin ice mantles. 



local ULIRGs (ISpoon et al 



20051 ). We do show some evidence for the 4.67 /im CO absorption 
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feature in three of these sources. Although tentative, this is the first suggestion of this feature 
in such high redshift systems. However this feature is very weak even in these sources and 
only seen in 3 of the 7 sources with usable spectra (and is not present in the stacked spectra). 
We can therefore conclude that it is either not ubiquitous even among sources with silicate 
absorption, or more realistically that the feature manifests with a range of absorption depths, 
and we only see the most pronounced features. 

As discussed in §4.5, a ~5/im starburst continuum would effectively weaken the ob- 
served CO feature depth (as that likely originates close to the AGN it self). For weak- PAH 
sources, the ~5/im continuum is most likely dominated by the AGN (ISajina et al.ll2007al ). 
however a starburst contribution, sufficient to affect the CO depth, cannot be excluded. 
The fact that one of our two strong-PAH sources (MIPS22432) is also one of three sources 
showing some hints of CO absorption, suggests that this is much more complicated that a 
simple correction scaled to the starburst strength (as revealed by the PAH). 

Drawing a parallel with the dust-embedded Galactic YSOs, perhaps there is a range from 
purely solid form CO to fully gaseous CO as the tempera tures rise and the CO previously 
deposited on grain surfaces sublimates (see IWhittetl 120031 . ; and references therein). Solid- 
phase CO shows a much narrower feature than its gas-phase counterpart, and in our low- 
resolution spectra would be undetectable, unless exceptionally strong. In local ULIRGs, this 
problem is present as well, and only the high resolution spec trum of IRAS08572+3915 has 
allowed to rule out CO ice in this source jGeballe et al-lbood ). To date, only NGC4945 has 
a CO ice detection (jSpoon et al.ll2003l ). 



Following the argument in iLutz et al.l (120041 ) , we might also conclude that we are not 
seeing this feature due to dust tori obscuration rather than 'buried nuclei'. However, this 
runs counter to our conclusions from the strong 3.4 /im HAC absorption features (see §5.3). 
Ultimately, there are too few 4.67 /xm CO detections even locally (see Introduction) and 
therefore the statistics on what fraction of the sources we should expect to show this feature 
sufficiently strongly are lacking. 



5.6. Buried nuclei at z ~ 2 



Our sources are drawn from sample of 24 /im flux limited and color-selected sources 
(jYan et al.l 120071 ) . The 0.9 mJy 24 /im flux limit already restricts us at 2~2 t o the highest 
lumin o sity sources, an d indeed we find (Lib) ~ 10 12 7 L for all our z ~ 2 sources (ISajina et al. 
20081 ). Ilmanishil (120091 ) shows that locally, at comparable luminosities, ~ 70% of the ULIRGs 
host buried AGN, compared to ~ 30% in lower luminosity ULIRGs. It is not surprising, 
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therefore, that in ISajina et al.l (12007a ) we found that the bulk of our z ~ 2 sources showed 
AGN-like mid-IR spectra, usually accompanied by silicate absorption features. Because of 
the biased selection of our sample, we cannot test whether or not obscured AGN form a larger 
fraction of z~2 ULIRGs than observed for the local ULIRGs. We did show that, consistent 
with the well-known luminosity evolution between z ~ 2 and today, the number density of 
our silicate-absorption sources is about an order of magnitude higher t han observed locally 
(for ULIRGs with comparable silicate absorption (ISajina et al.ll2007al )). However, silicate 
absorption alone is not a strong enough diagnostic of buried AGN since it can in principle 
arise from a range of different spatial scales (i.e. host obscuration). 

In this study, we find that beyond the presence of silicate absorption, the full 2-8 /im 
SEDs of our z ~ 2 sources, the presence and strength of 3.4 /im amorphous hydrocarbons and 
3.0 /im water ice absorption features are also comparable to those observed in nearby heavily 
obscured ULIRGs. This suggests similar obscuration geometries. Their closest counterparts 
in the local Universe are believed to be dominated by buried nuclei (i.e. compact sources with 
steep dust tempe rature gradients) mos t likely harboring AGN as discussed in the preceeding 



sections (see e.g. Ilmanishi et al.l 120071 . for a more complete discussion). The lack of strong 



4.67/zm CO absorption however is inconsistent with this picture, although the dearth of 
such observations locally makes it difficult to interpret (see §5.5). We find evidence of 
PAH emission even among the weak-PAH (i.e. AGN-dominated in the mid-IR) sources, 
suggesting that star-formation is present at least in some of these sources. Scaling from the 
3.3 /xm luminosities we expect that SFRs of order ~100M G /yr could be present. This is 
typically insufficient to dominate the bolometric output of these galaxies. The combination 
of deeply obscured, likely quasar, nuclei and ULIRG-strength star-formation supports the 
view of co-evolution of black-hole growth and host galaxy stellar mass growth suggested 



by a variety of observational and theoretical considerations (see e.g. iMagorrian et al.l 11998 



Canalizo fc Stockton 2001: Hopkins et al. 2008; Martmez-Sansigre et al. 2008). 



6. Summary & Conclusions 

1) In this paper we present deep rest-frame 2-8 /im spectra for 11 z~2 ULIRGs. 

2) Our sources tend to have redder 3-4/im continua than local ULIRGs. The few lo- 
cal sources with comparably red spectra, show continuum-dominated 2-10 /im spectra with 
weak-PAH and strong silicate absorption (consistent with their being predominantly AGN- 
dominated). They are therefore overall good analogs to our sources. 
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3) We found that nearly half of the sources with sufficient coverage (43%) show water 
ice absorption at 3/xm (compared with 65% for local ULIRGs). There is no clear trend of 
detectability with increasing silicate feature depth. 

4) The 3.4 HAC absorption is seen in 40% of the sample compared with 9% for local 
ULIRGs. The HAC-to-siilcate ratio, T3.4/V9.7, is 2-5 times stronger that the typical ratio 
observed in the Milky Way, but consistent with the obscured nuclei of some local ULIRGs. 
This is consistent with compact nuclei buried in dust with a steep temperature gradient 
distribution. 

5) Three of our sources show tentative signs of the CO 4.67 /xm feature. If genuine CO 
absorption, these features show a strong asymmetry between the short- and long-A branches 
of this feature. Moreover, the observed strength is still significantly lower than observed 
in local ULIRGs. Stacking analysis suggests this feature is not present in the bulk of our 
sample or is very weak (r 4 6 << 0.2). 

6) One source, MIPS22530, has a clear 3.3 /xm PAH emission feature detected. The 
ratio of the 6.2 to 3.3 /xm PAH strengths is consistent with that observed in local ULIRGs, 
and the overall strength of both features implies a starburst-dominated source, consistent 
with earlier results. Three more sources (MIPS 16080 MIPS22303, and MIPS22633) show 
tentative evidence for the 3.3 /xm PAH. The stacked spectra show evidence for PAH emission 
in the bulk of the weak- PAH sources implying likely significant star- formation in conjunction 
with the obscured black-hole accretion. 
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Table 1. Summary of observations 



MIPSID 


Zold 




G01 




G02 


G04 








SL1 


LL2 


SL1 


LL2 


SL1 


LL2 


464 


1.85 




6 x 120 s (S) b 






4 x 40 x 60s (M) 


30 x 120 s (S) 


8392 c 


1.9 


2 x 60s(S) 






12 x 120 s (S) 




30 x 120 s (S) 


15880 


1.68 




6 x 120 s (S) 






4 x 8 x 60s (M) 




16059 


2.325 




6 x 120 s (S) 






4 x 40 x 60s (M) 


30 x 120 s (S) 


16080 


2.007 




6 x 120 s (S) 






4 x 40 x 60s (M) 


30 x 120 s (S) 


16152 


1.83 








6 x 120 s (S) 


6 x 60s (S) 


30 x 120 s (S) 


22303 


2.34 




6 x 120 s (S) 






4 x 8 x 60s (M) 


30 x 120 s (S) 


22432 


1.59 






1 x 240 s(S) 


6 x 120 s(S) 






22530 


1.952 




6 x 120 s (S) 






4 x 40 x 60s (M) 


30 x 120 s (S) 


22548 


2.12 








6 x 120 s (S) 


6 x 60s (S) 


30 x 120 s (S) 


22633 


2.36 








6 x 120 s (S) 


4 x 8 x 60s (M) 


30 x 120 s (S) 



a Based on lSaiina et alJ l|2007al l and lDasvra et al.l j2009h ■ 

b We use the convention (^positions X )#cycles X ramp duration. In parentheses we mark whether the data is taken in 
stare mode (S) or spectral mapping (M). For simplicity the number of positions (i.e. the 2 nod positions) is omitted for 
stare mode observations. 

c This source is identified as FLS CV Z VLA9 from a GOl GTO program (PID#15, PI J. Houck). It is also identical to 
AMS06 llMartmez-Sansiere et al]|2008t) . 



Table 2. Continuum parameters 



MIPSID 






3-4 /im slope 


5 — 7/xm slope 


464 


11.07 


11.66 


3.12 


2.85 


8392 


11.63 


12.11 


3.63 


0.91 


15880 


11.31 


11.89 


2.52 


2.90 


16059 


11.74 


12.12 


3.15 


0.97 


16080 


11.49 


11.95 


3.38 


0.95 


16152 


11.68 


11.90 


2.12 


1.59 


22303 


11.95 


12.42 


3.30 


1.08 


22432 


11.41 


11.90 


3.75 


2.76 


22530 


11.10 


11.51 


1.90 


4.76 


22548 


11.88 


12.03 


2.26 


1.54 


22633 


11.54 


12.06 


3.37 


-0.30 
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Table 3. Optical depth of absorption features 



MIPSID 


Zold 


^ncw 


T3.05 


T3.4 


T9.7 


464 


1.85 


1.75 




0.52±0.10 a 


3.5±l.l b 


8392 


1.90 


1.90 




<0.45 


>4.9 


15880 


1.68 


1.62 




<0.28 


2.8 ±0.3 


16059 


2.325(Keck) 


2.325 


<0.22 


0.43 ±0.04 


1.9 ±0.6 


16080 


2.007(Keck) 


2.007 


0.33 ±0.13 


0.25 ±0.04 


1.5 ±0.4 


16152 


1.83 


1.84 


<0.28 


<0.16 


1.4 ±0.9 


22303 


2.34 


2.28 


0.80 ±0.19 


<0.24 


2.1 ±0.4 


22432 


1.59 


1.60 






1.8 ±0.6 


22530 


1.952 (Keck) 


1.952 


<0.34 




>3.7 


22548 


>2.12 


2.15 


<0.26 


<0.16 


1.5 ±0.8 


22633 


2.36 


2.20 


0.44 ±0.25 


0.37 ±0.14 


1.1 ±0.9 



a The errors on the ice and HAC features only reflect the per pixel errors in 
the spectra, assuming that errors in the continuum levels are negligible. We 
believe this is a reasonable assumption, in all cases except MIPS22633, as can 
be seen in Figure[2] 

b The erro r s on Tsi are based on the MCMC fitting as described in 
ISaiina et ail ll2007al) . They therefore include the per pixel errors, the un- 
certainty in the PAH+continuum decomposition, and the frequently poorly 
determined continuum above ~ 11 (im. For strong-PAH sources however, the 
depth of this feature is less well defined than shown here (see § 14.91 for further 
details). 



